Cosmology with tachyon field as dark energy 
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We present a detailed study of cosmological effects of homogeneous tachyon matter coexisting with 
non-relativistic matter and radiation, concentrating on the inverse square potential and the expo- 
nential potential for the tachyonic scalar field. A distinguishing feature of these models (compared 
to other cosmological models) is that the matter density parameter and the density parameter for 
tachyons remain comparable even in the matter dominated phase. For the exponential potential, 
the solutions have an accelerating phase, followed by a phase with a{t) oc t^^^ as t ^ oo. This elimi- 
nates the future event horizon present in ACDM models and is an attractive feature from the string 
theory perspective. A comparison with supernova la data shows that for both the potentials there 
exists a range of models in which the universe undergoes an accelerated expansion at low redshifts 
and are also consistent with requirements of structure formation. They do require fine tuning of 
parameters but not any more than in the case of ACDM or quintessence models. 



I. MOTIVATION 

Observations suggest that our universe has entered a 
phase of accelerated expansion in the recent past. Fried- 
mann equations can be consistent with such an accel- 
erated expansion only if the universe is populated by 
a medium with negative pressure. One of the possible 
sources which could provide such a negative pressure will 
be a scalar field with either of the following two types of 
Lagrangians: 
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Both these Lagrangians involve one arbitrary function 
V{(j)). The first one iquin, which is a natural general- 
ization of the Lagrangian for a nonrelativistic particle, 
L = (l/2)(7^ — V{q), is usually called quintessence (for a 
sample of models, see Ref. [1]). When it acts as a source 
in Friedmann universe, it is characterized by a time de- 
pendent w{t) = {PI p) with 
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Just as La 
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^qum generalizes the Lagrangian for the non- 
relativistic particle, itach generalizes the Lagrangian for 
the relativistic particle [2]. A relativistic particle with a 
(one dimensional) position q{t) and mass m is described 

It has the energy 
which 
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by the Lagrangian L = —ni\/l — cp- . 



and momentum p — mqj \f\ 



E = m/^/T- g2 

are related by = p'^+m''. As is well known, this allows 
the possibility of having massless particles with finite en- 
ergy for which E^ = p^. This is achieved by taking the 
limit of m ^ and g — + 1, while keeping the ratio in 
E = m/ y^l — (p- finite. The momentum acquires a life 
of its own, unconnected with the velocity g, and the en- 
ergy is expressed in terms of the momentum (rather than 



in terms of q) in the Hamiltonian formulation. We can 
now construct a field theory by upgrading q{£) to a field 
0. Relativistic invariance now requires (f> to depend on 
both space and time [0 — (j){t, x)] and q^ to be replaced 
by di4)d''4). It is also possible now to treat the mass pa- 
rameter m as a function of 0, say, V{(j)) thereby obtain- 
ing a field theoretic Lagrangian L = ~V{(I)) \J\ — d^4>9i4>- 
The Hamiltonian structure of this theory is algebraically 
very similar to the special relativistic example we started 
with. In particular, the theory allows solutions in which 

— > 0, di(j)d^<j) 1 simultaneously, keeping the energy 
(density) finite. Such solutions will have finite momen- 
tum density (analogous to a massless particle with finite 
momentum p) and energy density. Since the solutions can 
now depend on both space and time (unlike the special 
relativistic example in which q depended only on time), 
the momentum density can be an arbitrary function of 
the spatial coordinate. This form of scalar field arises 
in string theories [3] and — for technical reasons — is 
called a tachyonic scalar field. (The structure of this La- 
grangian is similar to those analyzed in a wide class of 
models called K-essence; see for example, Ref. [4].) This 
provides a rich gamut of possibilities in the context of 
cosmology [2,4-7]. 

The stress tensor for the tachyonic scalar field can be 
written in a perfect fluid form 



= ip+pyuk-p6i 



with 
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p^-v{(t>)^/T~d^^ 



The remarkable feature of this stress tensor is that it 
could be considered as the sum of a pressure less dust 
component and a cosmological constant [2]. To show this 
explicitly, we break up the density p and the pressure p 
and write them in a more suggestive form as 



P^Pv + Pdm; P = Pv+Pum, 



(5) 
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where 



V{(t))d'(t)d,(j) 



(6) 



PV = y(</>)^/^^ 



Pv = -pv 



This means that the stress tensor can be thought of 
as made up of two components - one behaving Uke a 
pressure-less fluid, while the other having a negative 
pressure. In the cosmological context, when (j) is small 
(compared to V in the case of quintessence or com- 
pared to unity in the case of tachyonic field), both these 
sources have w ^ —\ and mimic a cosmological con- 
stant. When (j)^ V ^ the quintessence has w k 1 leading 
to pg oc (1 + z)^; the tachyonic field, on the other hand, 
has w w for ^ ^ 1 and behaves like non-relativistic 
matter. 

An additional motivation for studying models based on 
-^quin is the following: The standard explanation of the 
current cosmological observations will require two com- 
ponents of dark matter: (a) The first one is a dust com- 
ponent with the equation of state p = contributing 
n„ 0.35. This component clusters gravitationally at 
small scales (Z ^ 500 Mpc, say) and will be able to ex- 
plain observations from galactic to super-cluster scales, 
(b) The second one is a negative pressure component with 
equation of state like p = wp with —1 < w < —0.5 con- 
tributing about fly ~ 0.65. There is some leeway in the 
{p/p) of the second component but it is certain that p 
is negative and (p/p) is of order unity. The cosmolog- 
ical constant will provide w = —1 while several other 
candidates based on scalar fields with potentials [1] will 
provide different values for w in the acceptable range. 
By and large, component (b) is noticed only in the large 
scale expansion and it does not cluster gravitationally to 
a significant extent. Neither of the components (a) and 
(b) has laboratory evidence for their existence directly or 
indirectly. In this sense, cosmology requires invoking the 
tooth fairy twice to explain the current observations. It 
was suggested recently [2] that one may be able to ex- 
plain the observations at all scales using a single scalar 
field with a particular form of Lagrangian. 

In this paper we explore the cosmological scenario in 
greater detail, concentrating on the background cosmol- 
ogy. Our approach in this paper will be based on the 
above view point and we shall treat the form of the La- 
grangian for Ltach as our starting point without worry- 
ing about its origin. In particular, we shall not make 
any attempt to connect up the form of V{(t)) with string 
theoretic models but will explore different possibilities, 
guided essentially by their cosmological viability. 

We construct cosmological models with homogeneous 
tachyon matter, assuming that tachyon matter co-exists 
with normal non-relativistic matter and radiation. Sec- 
tion II presents the equations for the evolution of scale 
factor and the tachyon field. In Section III we present 



solutions of these equations and discuss variations intro- 
duced by the available parameters in the model. Here 
we have analyzed two different models of the scalar field 
potential, one is the exponential potential and the other 
is the the inverse square potential which leads to power 
law cosmology [7]. In Section IV we compare the model 
with observations and constrain the parameters. Section 
V discusses structure formation in tachyon models. The 
results are summarized in concluding Section VI. 



II. TACHYON-MATTER COSMOLOGY 

For a spatially flat universe, the Friedman equations 
are 
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— P, - = - — (P + 3,) (7) 



where p = Pnr + Pr + with respective terms denoting 
non-relativistic, relativistic and tachyon matter densities. 
For the tachyon field we have 



P4' = 



y(0) 



p<t> 



(8) 



The equation of state for tachyon matter is p = wp with 



w 



1. The scalar field equation of motion is 



= -(1 - 



3H<j) + 
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V{cb) deb 



(9) 



The structure of this equation suggests that the change 
in goes to zero as it approaches ±1. In this case, the 
equation of state for the tachyon field is dustlikc. Thus 
at any stage if the tachyon field behaves like dust, it 
will continue to do so for a long time. This behavior 
persists for a duration that depends on the closeness of 
(j) to ±1. Detailed behavior will depend on the form of 
the potential. 

We shall discuss cosmological models with two different 
y(0) in this paper. The first one has the form 

''(« = 3^('-|;)"'*-'- 

It was shown in Ref. [7] that the above potential leads 
to the expansion a(t) = t" if is the only source. (The 
form V{<j)) DC 1/(0 — 0o)^ can be reduced to the form 
given above by a simple redefinition of the scalar field.) 
In this case, the term in the square bracket in equation 

(9) is 3H(j) — 2/(f) . This term vanishes in the asymp- 
totic limit when tachyons dominate and 4> rx t. This 
asymptotic solution is stable in the sense that all ini- 
tial conditions eventually lead to this state. If normal 
matter or radiation dominates, (j) stays close to the tran- 
sition point 2/{3H(j)), unless we start the field very close 
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to = 1. Thus the presence of non-rclativistic matter 
and radiation leads to a change in the equation of state 
for the tachyon field. The change in equation of state 
for the tachyon field implies that it is not a tracker field 
(for details see [8,9]). To get a viable model, i.e., matter 
domination at high redshifts and an accelerating phase 
at low redshifts, we need to start the tachyon field such 
that (j) is very close to unity, and (f) is very large. We will 
describe the fine tuning required in greater detail when 
we discuss numerical solutions of these equations. 

The second form of potential which we will consider is 
the exponential one with V{(f)) cx e-"^/"^". Then the term 
in square bracket in the equation (9) can be written as 

3H(p — 1/00 • In a universe dominated by radiation or 

normal matter, H{t) is a monotonically decreasing func- 
tion of time while (pQ is a constant. From this equation, 
it is clear that (p will increase if < l/{3H{t)(l)o). As 
H{t) keeps decreasing, will increase slowly and ap- 
proach unity, asymptotically. However, in the mean- 
while, tachyons may begin to dominate in terms of energy 
density and this changes the behavior. In a tachyon dom- 
inated scenario, we get rapid expansion of the universe 
and H{t) varies much more slowly than in the matter or 
radiation dominated era. Thus </> changes at a slower rate 
but it still approaches unity and hence we — eventually 
— get a dust like equation of state for the tachyon field. 
Whether there is an accelerating phase or not depends on 
the initial values of (p, <po and fl^. Present values of den- 
sity parameter for non-relativistic matter and tachyons 
fixes the epoch at which tachyons start to dominate the 
energy density. The parameter (Pq sets the time when (p 
approaches unity and the asymptotic dust like phase for 
tachyons is reached. Initial value of </> fixes the duration 
of the accelerating phase. If this value is very close to 
unity then it departs very little from this value through 
the entire evolution. On the other hand, if it starts far 
away from this value, then the equation of state for the 
tachyon field can lead to a significantly long accelerat- 
ing phase. It is possible to fine tune the evolution by 
choosing a sufficiently large value for the constant (po, so 
that (p is much smaller than unity even at the present 
epoch, and by requiring that the tachyon field is start- 
ing to dominate the energy density of the universe at the 
present epoch. In such a case, the accelerating phase is 
a transient between the matter dominated era and the 
tachyon dominated era with a{t) cx i^/^ in both of these 
regimes. This model has been the focus of many studies 
[10] though most of these have chosen to ignore the role 
of matter or radiation. 

Some of the models discussed below have an acceler- 
ating phase for the universe followed by a dust-like ex- 
pansion with a(t) cx t^/^ asymptotically. These mod- 
els have the attractive feature that they do not be- 
have as a do Sitter-like universe in the asymptotic future 
and thus do not possess a future horizon. String theo- 
retic models have difficulty in incorporating cosmologies 
which approach de Sitter-like phase asymptotically; from 



this point of view, these models can be theoretically 
attractive. 



III. NUMERICAL SOLUTIONS TO 
COSMOLOGICAL EQUATIONS 

This section presents a detailed study of the back- 
ground tachyon cosmology. We solve the cosmological 
equations of motion numerically for the two potentials 
mentioned in the previous section. 



A. Inverse square potential 

We start with the discussion of the potential V{(p) cx 
Energy density of the universe is a mix of the 
tachyon field, radiation and non-relativistic matter. The 
asymptotic solutions for the universe with only the 
tachyon field as source are well understood: we get rapid 
expansion of the universe with the field (p growing in pro- 
portion with time. We wish to see the effect of other 
species (matter and radiation) on the evolution of the 
tachyon field, and the stability of the asymptotic solu- 
tion. 

Since only two of the three equations in (7) and (9) are 
independent, we choose to drop the equation with second 
derivative of the scale factor. We choose some instant of 
time as tin and rescale variables as follows, for numerical 

convenience: 



tHi. 



y = a{t)/a{ti„), y' = dy/dx 



(10) 



and ip = (p/(pin, ip' = dtp/dx. (We will use the present 
epoch and Planck epoch for tin, in two different contexts 
described below.) Then the equations are 



y =y 
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The subscript in refers to the initial value, and f^M and 
fifl refer to the density parameters for matter and ra- 
diation respectively. It is assumed that the sum of all 
density parameters is unity. 

The initial conditions for y and follow from their 
definition: yin = 1, ipin — 1- The initial condition for 
tp', can be related to the density parameter of tachyon 
matter at t = tin 
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FIG. 1. Phase portrait for the scale factor. This plot is for 
f^M (present) = 0.3 and n = 6. We choose the present value of 
4>inHiri, — 2.4 where 'in' refers to present epoch [a — 1; z — 0) . 
The transition from a decelerated expansion to an accelerated 
one is clear in this picture. We have used this model for all 
the figures pertaining to the inverse square potential. 

There are two branches of solutions, one with positive 
tp' and another with negative values of For each set 
of ^mi„ and n, there is a minimum value of (j)inHin below 
which wc cannot satisfy this equation. We can solve these 
equations if tp'^ > 0, or if 
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3n 
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For a given value of (jxf^Hf^, we get initial values of all 
the variables from the above equations and using these 
wc evolve these quantities. As shown in above equations, 
parameters n, il^^^, 4>inHin and \(j>in\ are inter-related. 
For a fixed n and 0^^^, larger values of (pf^Hf^ imply a 
larger value of initial |(/)|, and hence an equation of state 
that is closer to that of dust. 

We study the evolution of such a model in three steps. 
First wc evolve the system from present day to future and 
show that the asymptotic solution [7] is stable, i.e., the 
entire allowed range of present values leads to the same 
asymptotic solution. In this context 'in' refers to the 
present values of the parameters. We choose the present 
value of r^Min = 0.3, and all the results shown here use 
n — 6 (results do not change qualitatively for other values 
of n), i.e., the asymptotic solution gives a{t) oc t^. 

Wc plot the phase portrait for the scale factor in Fig. 1. 
This figure clearly shows the late time acceleration in 
tachyon dominated universe. (The figures exhibit past as 
well as future evolution; the past evolution is discussed 
later.) The positive branch is plotted in this figure and 
the following ones but the discussion holds true for the 



negative branch as well. Evolution of field cf) and its time 
derivative (f> is shown in Figs. 2 and 3. These too demon- 
strate the approach to the asymptotic solution for the 
two initial conditions. In asymptoti c futu re, we expect 
to approach a constant value of y^2/3n. Indeed, both 
the negative and the positive branches (which respec- 
tively correspond to positive ip' and negative ip') exhibit 
this behavior. The fact that other initial conditions lead 
to similar behavior is shown in the phase diagram for 
field (j) in Fig. 4. We see that for the entire range of ini- 
tial conditions for the field (p, it quickly approaches the 
asymptotic solution. It takes longer to reach the asymp- 
totic solution for initial conditions where \(j)inHin\ (the 
value at present) is much larger than its minimum al- 
lowed value as in this case 101 is closer to unity and hence 
the present equation of state for the tachyon field is like 
that of a pressure-less fluid. In such a case accelerating 
phase of the universe starts at late times compared to 
models where \(l)inHin\ is closer to its minimum allowed 
value. Thus we need to fine tune the value of this pa- 
rameter in order to arrange for accelerating phase of the 
universe to start at low redshifts, as required by obser- 
vations. 




l+z 

FIG. 2. The tachyon field 0/0present as a function of red- 
shift. The model is the same as in Fig. 1 

In future evolution, density parameter for matter de- 
creases with the start of the accelerating phase. The 
density parameter for tachyons begins to increase in this 
phase and quickly approaches unity. 

Next we study the behavior of these solutions in past, 
i.e., we use the present day conditions and evolve back 
to sec the kind of initial conditions that arc required for 
us to get a viable model today. Figures used to illustrate 
future evolution also show the past evolution of these 
quantities. As motivated in section II, wc expect mat- 
ter and radiation to drive <p away from unity towards 
smaller values in forward evolution. So, as we evolve 
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the equations back, we expect all present day conditions 
to lead towards 4> = 1. This is indeed what we find 
in numerical solutions. The equation of state for tachyon 
field approaches that of a pressure-less fluid, and the uni- 
verse expands with a{t) oc t"^^^ at high redshifts. At even 
earlier times, radiation takes over and both matter and 
tachyons become subdominant components. 
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FIG. 3. This plot shows as a function of redshift for the 
same model. 




^/ ^Present 

FIG. 4. Phase plot for the tachyon field. Here, we started 
from initial conditions with a fixed value of and n. (p was 
varied and fi^.^ was allowed to vary with it. Both the positive 
branch as well as the negative branch are shown. The scale 
factor a{t) = 1/(1 + 2:) is marked along the arrow heads. 
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FIG. 5. SIm, and VLr as functions of redshift. The mat- 
ter density parameter is almost a constant at 2 < « < 10^ and 
then drops to small values as radiation begins to dominate. 
The present day values of density parameter for matter and 
radiation are 0.3 and 10~® respectively. 

The time dependent matter density parameter 
^^Af(i) = 8iTGpM{t)/-iH'^{t) is plotted as a function of 
redshift in Fig. 5. Here non-relativistic matter dom- 
inates in the sense that Qm > ^4, but evolution does 
not drive fl^ towards zero or Hm close to unity. Instead 
the two density parameters remain comparable during 
much of the "matter dominated" era (from z w 10^ to 
z « 2). This is a unique feature not seen in other cosmo- 
logical models. In the phase where the equation of state 
of tachyon field approaches w ~ 0, the ratio of density 
parameter for non-relativistic matter and tachyon field 
becomes a constant. For models of interest, i.e., where 
accelerating phase is starting at low redshifts, this ra- 
tio is of order unity for most choices of parameters and 
the relative importance of tachyons and matter does not 
change up to the time when accelerating phase begins. 
(This has a significant effect on gravitational instability 
and growth of perturbations; see section V.) 

The tachyon field (f> approaches a nearly constant value 
in past and changes very little from its value in the early 
universe to the time when the accelerating phase begins. 
In the pre-acceleration phase, <p is close to unity. This 
means that the field (f> starts from large values in the 
models that satisfy basic observational constraints. This 
is the fine tuning required in constructing viable tachyon 
models. In Fig. 6 we show the evolution of 4'{t)H{t) as a 
bmction of redshift. As expected, this is a constant in the 
asymptotic regime where tachyons dominate. However 
this function evolves rapidly in the dust like regime. As 
<f> remains a constant in this regime, the change is mainly 
due to change in H{t). The asymptotic value of this 
product is 4>H(t) — > -\/2n73, so we expect (t>H{t) to 
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decrease from its initial value and make a transition to 
the constant asymptotic value. We need to choose the 
initial value of this product so that the transition from 
dust like to accelerating phase happen s arou nd now. As 
the asymptotic value is (j)H{t) — > A/2n/3, and this is 
of order unity for small n, the value of 4> in the early 
universe has to be large. The product (j){tpi)H{tpi) is 
approximately (f>{to)H{tpi) as (j) does not change by much 
during its evolution. Since the present day (j){tQ)HQ is 
of order unity, (j){tpi)Hpi « Hpi/Hq. (Here to is the 
present time and Hq is the present value of the Hubble 
constant.) Thus the initial value of (j> sets the time scale 
when transition to asymptotic behavior takes place. This 
fine tuning is similar to that needed in models with the 
cosmological constant as the value of A has to be tuned 
to arrange for accelerating phase to start at low redshifts. 




0.1 1 10 100 1000 10* 10^ lo' 
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FIG. 6. Plot of (j)H as a function of redshift. 

To illustrate this point further, we start the integra- 
tion at Planck time for a range of initial conditions and 
show that indeed we obtain this behavior. The age of 
the universe at the time of transition from pressure-less 
behaviour for the tachyon field to an effective cosmolog- 
ical constant behavior is indeed proportional to the ini- 
tial value of (p. To illustrate this, we first plot w for the 
tachyon field as a function of scale factor in Fig. 7. This 
is plotted for many values of (t){tpianck)H{tpianck)- An or- 
der of magnitude increase in the initial value delays the 
onset of accelerating phase by roughly an order. Hence to 
achieve accelerating phase at late times, one requires the 
fine tuned value mentioned above. The same behavior 
is summarized in Fig. 8, where we plot the scale factor 
when = w^^^^^^^^^^^jl = (2/(3n) — l)/2 as a function 

of (t){tplanck)H{tpianck)- 




a(t)/a(t_Pl) 



FIG. 7. Plot of = (A^ — 1 as a function of the scale factor. 
The curve on the left has the lowest value of initial (pH and 
increases by an order of magnitude for respective curves as 
we go towards right. 




a(t)/a(tpi) 

FIG. 8. Plot of 4>inHi„ (at Planck time) as a function of 
the scale factor value at which = — 1 becomes half its 
asymptotic value. 



B. Exponential Potential 

We now repeat the above analysis for the exponen- 
tial potential. As before, we consider a mix of non- 
relativistic matter, radiation/relativistic matter and the 
tachyon field. We fix Q.^ = 0.7 at the present epoch. (In 
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the literature, it has been suggested that the exponential 
potential does arise in some of the string theoretic mod- 
els. However, in our approach, we think of V{4>) as an 
arbitrary function just as in quintessence models.) For 
the purpose of numerical work, we have taken the value 
of (f)/(j)o to be 10^ at the present epoch. The choice of (j) 
fixes the amplitude of the potential as we have already 
fixed flff, and <p/(f)o. The remaining parameter is 0o, and 
as we have already fixed 4'in/4>a, fixing 0o is equivalent to 
fixing (j)in. We use the combination (pinHin to construct a 
dimensionless parameter which indicates the value of (/)o. 
4'inHin should bc of Order unity if (f>in is comparable to 
the age of the Universe today. We shall study variation 
of quantities of interest with time for a range of values of 
<j) and 4'inHin. 

In the case of exponential potential, the equation of 
state for the tachyon field approaches that of dust as 
f — > cxD. Thus any accelerating phase that we might get 
will be followed by a dust like phase. The accelerating 
phase occurs when the tachyon field dominates the energy 
density of the universe as it is only in this case that we 
can get /> + 3p < 0. For an accelerating phase to exist 
at all, the tachyon field should begin to dominate before 
it enters the dust like phase. Thus the duration of the 
accelerating phase will depend on how different 4> is from 
unity at the start of tachyon dominated phase. 
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FIG. 9. Phase plot for the exponential potential. The du- 
ration of the accelerated phase depends on the initial value 
of 0. The accelerated phase ends into a dust like phase 
of tachyons. Solid lino is for = 0.1, dashed lino is for 
4>in = 0.3, dot-dashed, dashed and dot-dot-dashed curves are 
for 0i„ = 0.5, 0.7 and 0.9 respectively. 
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FIG. 10. Plot of density parameters for matter and tachyon 
as functions of redshift. The models are the same as in the 
previous figure. 

The phase plot for the scale factor is shown in Fig. 9. 
We have kept ^inHin fixed at present day and varied ^ 
for a range of values. The unique feature of this model 
is that we regain the dust like phase in the future. It 
is clear from this plot that the duration of accelerating 
phase varies considerably across models shown here. The 
supernova observations require the universe to be in an 
accelerating phase at z < 0.25, so models that do not 
have an accelerating phase at all or too small an acceler- 
ating phase can be ruled out easily. 

The fact that these cosmological models have two de- 
celerating phases, with an accelerating phase sandwiched 
in between, is noteworthy. In such models, one can ac- 
commodate the current acceleration of the universe with- 
out the model "getting stuck" in the accelerating phase 
for eternity. Since the universe has a{t) oc f^/'"* asymptot- 
ically, it follows that there will be no future horizon (for 
some other attempts to eliminate the future horizon, see 
[11]). String theoretic models have difficulty in accom- 
modating an asymptotically de Sitter-like universe and 
our model could help in this context. 

The plot in Fig. 10 shows Qm, ^tt> as functions of red- 
shift and shows that f2^ increases monotonically with 
time. 

A plot of as a function of fi^ in Fig. 11 shows that 
all initial conditions lead to = 1. The function (j)H 
increases with 1 + z. It is clear that all initial condi- 
tions lead to a narrow range of asymptotic values for this 
product. 
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FIG. 11. Plot of (j) Vs. $7^. The asymptotic value = 1 is 
reached irrespective of the initial conditions. The models are 
the same as in the previous figure. 

If we keep the value of (f) fixed and increase (pinH^n, 
the duration of the accelerating phase increases. This is 
illustrated in Fig. 12. A plot of the density parameters 
is shown in Fig. 13. 




0.01 0.1 1 



a(t) 

FIG. 12. Phase plot for a fixed value oi (j) = 0.38 and vary- 
ing (pinHin (present day). The solid lino is for <j)inHin = 1/3, 
dashed line, dot-dashed line, dotted and dot-dot-dashed line 
for (j)inHin = 5/3, 7/3, 3 respectively. As we increase the value 
of (f>inHin the duration of accelerated phase increases. 



1+z 

FIG. 13. This plot shows density parameters for matter 
and tachyon for a fixed value of 4> = 0.38 and varying (pinHin- 
The models are same as in Fig. 12. 

IV. COMPARISON WITH SUPERNOVA lA 
OBSERVATIONS 

Wc now compare the results with supernova la obser- 
vations by Perlmutter et al. [12,13]. The data from the 
redshift and luminosity distance observations of these su- 
pernovae, collected by the Supernova Cosmology Project, 
concludes that the universe is currently going through 
an accelerating phase of expansion. This is a powerful 
constraint on cosmological models, e.g. it rules out the 
Einstoin-deSitter model at a high confidence level. The 
obvious contender is a positive cosmological constant, 
though more exotic matter driving the accelerated ex- 
pansion is also a possibility. Since currently the energy 
density of the universe is dominated by a form of matter 
with negative pressure, the presence of a tachyon source 
is thus expected to be consistent with these observations. 

The results of theoretical model and observations are 
compared for luminosity measured in logarithmic units, 
i.e., magnitudes defined by 

meiz) = M + Uogw{DL) (15) 

whereAl = M — 5logio{Ho) and Dl = HqcIl, the factor 
M being the absolute magnitude of the object and dL is 
the luminosity distance 

dr. = {l + z)a{to)r{z); riz) ^ c J ^ (16) 

where z is the redshift. 

We perform the test of goodness-of-fit on the model 
with V{^) oc For our analysis, we have used all 



the 60 supernovac quoted in [12]. Of these 42 are high 
red shift supcrnovae reported by the supernova cosmol- 
ogy project [12,13] and 18 low red shift supernovae of 
Callan-ToUolo survey [14,15]. We have three parameters: 
flM{in), n and (t)inHin (with 'in' referring to present 
epoch). An additional freedom is the choice of sign of 
(j){x = tinHin). We freeze ^luiin) at 0.3 and we find that 
that for any choice of n, we can get a reasonably low value 
for reduced by choosing 4>oHq judiciously. In addition, 
the value of does not vary much over the range of pa- 
rameter values that we have studied. Minimum value of 
X^ per degree of freedom that we encountered is aroimd 
1.93. This test does not isolate any particular region in 
the parameter space so we shall refrain from quoting any 
particular values of parameters as our best fit. Suffice 
to say that a large range of each of these parameters is 
allowed by the supernova observations. This discussion 
holds for both the potentials discussed above. 

We have plotted distance modulus dm{z) = m — M as 
a function of redshift for one of the models in Fig. 14. 
The data points for the 60 supernovae are over plotted 
as well. The contours of reduced x^ and the value of 
ilm at redshift z = 10 are illustrated in Fig. 15. This 
is of interest while studying structure formation in these 
models. 
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6. H. /<i. ,H. , 

FIG. 15. The figure shows contours of Qm at z = 10 
(dashed lines) and contours of reduced (solid lines) for 
V{cl)) oc 1/0^. There is small region where a small 
overlaps Qm > 0.9. This region is near n = 1.2 and 
\og{(j)inHin/4>iniHini) = 0.05 The values along the x axis 
arc ratios of <j}inHi„ to the minimum value of (j>inHi„ for a 
particular value of power n. 



1 + z 

FIG. 14. Comparison of the model with n 
present (j)H = 2.56 with supernova la data. 
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We also did a x^ analysis of models with a range of 
parameters for the exponential potential. Here again the 
theoretical models satisfy the supernova constraints and 
we cannot rule out a specific model by this analysis. We 
plot the contours oOIm at z = 10 and x^ in Fig. 16. We 
obtain results similar to those for 1/0^ potential. 



FIG. 16. Contours of Qm at « = 10 (dashed lines) and 
contours of reduced (solid lines). The favored region is 
1 and 0.25 < d) < 0.4. 
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V. STRUCTURE FORMATION IN TACHYONIC 
MODELS 

Cosmological models with tachyons 

and non-relativistic matter have a significantly different 
behavior as compared to quintessence or the cosmologi- 
cal constant models. The most important difference here 
is that the source of acceleration in the Universe makes 
an insignificant contribution to the energy density of the 
universe beyond z > 1 in quintessence or cosmological 
constant models, whereas in tachyon models the density 
parameter for the scalar field docs not become insignifi- 
cant in comparison with the density parameter for mat- 
ter in most models. This has important implications for 
structure formation in these models as the density pa- 
rameter of the matter that clusters is always smaller than 
unity, and the rate at which perturbations grow will be 
smaller than in standard models. The exception to this 
rule is a small subset of models where fi^ approaches val- 
ues much smaller than unity beyond z ^ 1. As can be 
seen in the contour plots in the previous section, these is 
a small subset of the models that satisfy the constraints 
set by supernova observations. 

Given that the density parameter for matter is almost 
a constant at high redshifts {3 < z < 10^), we can solve 
for the rate of growth for density contrast in the linear 
limit. The equation for the density contrast is given by 
[16] 



2-S = AirGpS 

a 



(17) 



where 6 = {p — p)/p, the factor p being the average den- 
sity. Rescaling in the same manner as the cosmological 
equations we have 



(18) 



Here flMo is the density parameter for non-relativistic 
matter at the epoch when the scale factor is ao and 
the Hubble parameter is Hq. Since, Qm is nearly 
constant at high redshift, we get 5 oc t"^, and m = 
(1/6) (Vl + 24nM — 1) for the growing mode. The unique 
feature here is that at high redshift, matter density pa- 
rameter does not saturate at unity for all the models. 
This is true for both exponential and 1/0^ potentials. 
For models where the matter density parameter does not 
reach unity, the growth of perturbations is slow, as can 
be seen from the above equation. The models in which 
density parameter is unity at high redshifts, the growth 
of perturbations is closer to that in the ACDM model. 
The slower growth of perturbation implies that rms fluc- 
tuations in mass distribution were larger at the time of 
recombination as compared to conventional models. This 
will have an impact on the temperature anisotropies in 
the microwave background in these models. Since the 
latter is tightly constrained from CMBR measurements, 
models with slow growth of perturbations can be ruled 
out. 
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FIG. 17. Evolution of the density contrast with redshift. 
The curves correspond to f2i„ = 0.3 (present day value). The 
solid curve is for ACDM model given for comparison. The 
density parameter Om « 0.98 at a = 10 for the dashed line 
and Q,M « 0.56 at ^; = 10 for the dotted line. 

We use the solution outlined above to set the initial 
conditions for equation (18) and evolve forward through 
the regime where tachyons begin to dominate and mat- 
ter becomes irrelevant. As the universe begins to accel- 
erate, at late times wc anticipate that the growing mode 
would slow down and eventually saturate. This is in- 
deed true. The rate of growth for 6 slows down once the 
universe begins to accelerate. It comes to a halt around 
the epoch where the accelerating phase begins. This late 
time behavior is similar to what happens in most mod- 
els where the universe begins to accelerate at late times. 
The evolution is illustrated in Fig. 17. Here we have plot- 
ted two different models, one in which density parameter 
saturates at a small value and the other in which it ap- 
proaches unity. 

It is clear from the figure that one can indeed construct 
models in which the growth of perturbation is very sim- 
ilar to that in ACDM models. Such models are clearly 
viable. It is also obvious that these models are confined 
to a narrow range of parameters, as described in the fig- 
ure captions; if one moves out of this range, then the 
perturbations grow more slowly and should have higher 
amplitude in the past in order to maintain a given ampli- 
tude today. These are ruled out by CMBR observations. 
The following caveat, however, needs to be kept in mind 
in ruling out such models. It was suggested in [2] that one 
can construct models with tachyonic scalar field in which 
the equation of state is different at small scales and large 
scales. In such models, our conclusions will apply only at 
large scales and growth of structure at small scales will 
still be possible, i.e., inhomogeneities in the tachyon filed 
will play and role and may offset the conclusions about 
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growth of perturbations at small scales whereas our re- 
sults for the expansion of the universe will remain valid 
at sufficiently large scales. 

VI. CONCLUSIONS 

We have shown that it is possible to construct viable 
models with tachyons contributing significantly to the 
energy density of the universe. In these models, matter, 
radiation and tachyons co-exist. We show that a subset 
of these models satisfy the constraints on the accelerating 
expansion of the universe. For the accelerating phase to 
occur at the present epoch, it is necessary to fine tune 
the initial conditions. 

We have further demonstrated that the density param- 
eter for tachyons does not becomes negligible at high 
redshifts, hence the growth of perturbations in non- 
relativistic matter is slower for most models than, e.g., 
the ACDM model. This problem does not effect a small 
subset of models. However, given that the density pa- 
rameter of tachyons cannot be ignored in the "matter 
dominated era" , it is essential to study the fate of fluc- 
tuations in the tachyon field. 
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